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Exploring the QCD phase diagram with compact stars 
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We investigate a nonlocal chiral quark model with separable 4-fermion interaction for the case of (7(3) flavor 
symmetry and show that strange quark matter is unlikely to occur in a large enough volume of a compact star 
to entail remarkable observational consequences. The phase diagram in the two-flavor sector of such model has a 
critical end point of the line of first order chiral/deconfinement phase transitions on which a triple point marks 
the junction with the critical line for second order phase transitions to two-flavor color superconductivity (2SC) 
below T ~ 80 MeV. Stable hybrid star configurations with large quark matter core in a color superconducting 
phase can exist. A consistent cooling phenomenology requires that all quark species be gapped, the minimal 
pairing gap of the order of 10 — 100 keV. 



1. Introduction 

The investigation of the QCD phase diagram 
has become a research topic of highest prior- 
ity. Relativistic heavy ion collisions in the RHIC 
era have provided results which confirm previ- 
ous information about the critical temperature 
T c ~ 170 MeV of the hadronization transition in 
the approximately baryon-free regime pp, which 
parallels the situation a few microseconds after 
the big bang. This value is in agreement with the 
deconfinement transition temperature calculated 
in 2+1 flavor Lattice QCD These simula- 
tions are now extended into the region of (small) 
finite chemical potentials where one of the char- 
acteristic features is a critical endpoint of first 
order phase transitions (tricritical point). Future 
heavy-ion collision experiments planned at GSI 
Darmstadt (FAIR) will explore the phase diagram 
in the finite density domain and hope to find ex- 
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perimental signatures of the tricritical point. In 
dense quark matter at temperatures below ~ 50 
MeV, due to attractive interaction channels, the 
Cooper pairing instability is expected to occur 
which should lead to a variety of possible quark 
pair condensates corresponding to color supercon- 
ductivity (CSC) phases, see jS] for a review. Since 
it is difficult to provide low enough temperatures 
for CSC phases in heavy-ion collisions, only pre- 
cursor phenomena |4I5| are expected under these 
conditions. CSC phases may occur in neutron 
star interiors |Hj and could manifest themselves, 
e.g., in the cooling behavior |7I8| . 

However, the domain of the QCD phase dia- 
gram where neutron star conditions are met is 
not yet accessible to Lattice QCD studies and 
theoretical approaches have to rely on nonpertur- 
bative QCD modeling. The class of models clos- 
est to QCD are Dyson-Schwinger equation (DSE) 
approaches which have been extended recently to 
finite temperatures and densities Within sim- 
ple, infrared-dominant DSE models early studies 
of quark stars ^H] and diquark condensation JI] 
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have been performed. 

The present contribution will be based on re- 
sults which have been obtained within a nonlo- 
cal chiral quark model (NCQM) using a covari- 
ant, separable representation of the gluon propa- 
gator ^2] with formfactors suitable for an exten- 
sion to finite temperatures and chemical poten- 
tials |13ll4ll5j . Alternative approaches to sepa- 
rable chiral quark models are based on the instan- 
ton approach [16117118] . As a limiting case, the 
well-known Nambu-Jona-Lasinio (NJL) model is 
obtained for a cutoff formfactor. For recent appli- 
cations of the NJL model to QCD at high density 
in compact stars, see and references therein. 
A different class of QCD-models is of the bag- 
model type. These models do not describe the 
chiral phase transition in a self-consistent way 
and will not be considered here. We will con- 
sider within a NCQM the following questions: (i) 
Is strange quark matter relevant for structure and 
evolution of compact stars? (ii) Are stable hybrid 
stars with quark matter interior possible? (iii) 
What can we learn about possible CSC phases 
from neutron star cooling? 

2. Nonlocal chiral quark model, Nt = 3 

In order to answer the question whether strange 
quark matter phases should be expected in the 
neutron star interiors, we consider here a three- 
flavor generalization of a NCQM with the action 



S[q,q] 



+Dr 



K Q=0 



{q{k)(Jc - m)q(k) 



where J k — J ■ We restrict us here to the 
scalar current jf (k) — q(k)\ a f(k)q(k) and the 
pseudoscalar current jp{k) — q{k)i^X a f{k)q(k) 
in Dirac space with q(k) and q(k) being quark 4- 
spinors and the formfactor f(k) accounts for the 
nonlocality of the interaction. For the applica- 
tions to compact stars, the quark matter equation 
of state is needed and can be obtained from the 
partition function 



Z [ T , A] = J VqVq exp \ S[q, q] - fi^qq^j , 



(2) 



where the constraint of baryon number conserva- 
tion is realized by the diagonal matrix of chem- 
ical potentials jl (Lagrange multipliers). Fur- 
ther details concerning the notation and the 
parametrization of this model can be found in 
Ref. EOj. 

In order to perform the functional integrations 
over the quark fields q and q we use the formalism 
of bosonisation which is based on the Hubbard- 
Stratonovich transformation of the four-fermion 
interaction. The resulting transformed partition 
function in terms of bosonic variables will be con- 
sidered in the mean-field approximation 



Z[T,fi] 



/ "7 



f ' 



(3) 



with the 4- vector kf — (ko — /i/, k) and the effec- 
tive quark masses Mf — Mf(k) = to/ + <fiff(k) 
containing the flavor dependent mass gaps <pf. In 
the mean field approximation, the grand canoni- 
cal thermodynamical potential is 

tt(T,M) = yhx{Z[T,{» f }]/Z[0,{0}]} , (4) 

where fluctuations are neglected and the diver- 
gent vacuum contribution has been subtracted. 
The quark mass gaps 4>f are determined by solv- 
ing the gap equations which correspond to the 
minimization conditions dil/d(f>f = and read 



8D N C 



d A k 2M / /(fc) 

(27T) ^ 



h-2 



(5) 



As can be seen from ©, for the chiral t/(3) quark 
model the three gap equations for <j> u , <f>d, 4> s are 
decoupled and can be solved separately. 

In what follows we consider the case T = 
only. All thermodynamical quantities can now 
be derived from Eq. Q, for details see [20]. For 
instance, pressure, density and the chiral conden- 
sate are given by 



p = -SI , n } = -- — , <qq>= -- 



dfi f 



(6) 



the energy density follows from e = — p + 
X)// 1 / 71 /- For most of the numerical investi- 
gations within the NCQM described above we 



3 



will employ a simple Gaussian formfactor f(k) = 
exp(fc 2 /A 2 ) , which has been used previously for 
studies of meson and baryon properties in the 
vacuum |17| as well as quark deconfinement and 
mesons at finite temperature |13I15| . A system- 
atic extension to other choices of formfactors can 
be found in JSj- The Gaussian model has four 
free parameters to be defined: the coupling con- 
stant Do, the interaction range A, and the current 
quark masses m u = nid, rn s . These are fixed by 
the pion mass = 140 MeV, kaon mass iuk = 
494 MeV, pion decay constant f„ = 93 MeV and 
the chiral condensate — < qq > = (240 MeV) 4 . 
Fig. ^ shows the behavior of the thermodynami- 




Figure 1. Dependence of the thermodynamical 
potential on the light flavor gap <p q — <j) u — 4>d 
(order parameter) for different values of the chem- 
ical potential, <f) s — 682 MeV. 



cal potential as a function of the light quark mass 
gap <j) q — 4> u = 4>d for different values of the chem- 
ical potential /i. For /i < \i c — 330 MeV the argu- 
ment and the value of the global minimum is in- 
dependent of fj, which corresponds to a vanishing 
quark density (confinement). At /i = \x c a phase 
transition from the massive, confining phase to a 
deconfining phase with negligibly small mass gap 
occurs. In Fig. [5] we show that the quark pres- 
sure vanishes for /i < /i c (confinement) and can 
be well described by a two-flavor bag model with 
a bag constant B = 81.3 MeV/fm 3 in the region 
of chemical potentials 330 MeV < fx < 492 MeV 
where the third flavor is still confined. 



Separable model will] 3 flavors 

2 11. Bag model B= 81.3 MeV/fm 3 

3 11 Baa model B= 100.7 MeV/fm" 1 



Figure 2. Pressure of the NCQM as a function 
of the chemical potential for the separable model 
(solid line) compared to a three-flavor (dotted 
line) and a two-flavor (dashed line) bag model. 
All models have the same critical chemical po- 
tential fi c ~ 330 MeV for (light) quark deconfine- 
ment. 

3. Color superconductivity, Nf = 2 

The quark-quark interaction in the color anti- 
triplet scalar diquark channel is attractive driving 
the pairing with a large zero-temperature pairing 
gap A ~ 100 MeV for the quark chemical poten- 
tial /jL q ~ 300 — 500 MeV. Therefore, we consider 
now a NCQM described by the effective action 

3A) 



S[q,q] 



k.k' 



jl(k)j d (k'), (7) 



where jd(k) = q 7 \k)Ci"f^T2\2f{k)q{k) is the 
scalar diquark current with the charge conjuga- 
tion matrix C = ^7072- After bosonization, the 
mean field approximation introduces a new or- 
der parameter: the diquark gap A, which can be 
seen as the gain in energy due to diquark conden- 
sation. The mass gaps 4> u , <fid indicate dynamical 
chiral symmetry breaking. The grand canonical 
thermodynamic potential per volume can be ob- 
tained as 



a, 



,A: 



Vd, T) 



D 



2 ^ 



d 3 k 



-Tr In 



— G 1 {iuj„,k) 



3 D 

, (8) 



where io n = (2n + \)nT are the Matsubara fre- 
quencies for fcrmions and the inverse Nambu- 
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Gorkov quark propagator is 

r -i = ( ft- M - £7o A 75 t 2 A2/(£:) \ 
~ \ -A* 75 T 2 A 2 /(fc) ft - M + /no J 
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Figure 3. Phase diagram for Nf = 2 quark matter 
in the NCQM. 

The resulting phase diagram is shown in Fig. [21 It 
includes a 2-fiavor color superconductivity (2SC) 
phase for which quarks of one color, say blue, re- 
main unpaired. The color-flavor locking (CFL) 
phase HJ| requires approximate SU(3) flavor sym- 
metry and can be excluded from our discussion 
since strange quarks remain confined up to the 
highest densities occuring in a compact star con- 
figuration pu) . 

4. Structure and cooling of hybrid stars 

In describing the hadronic shell of a hybrid star 
we exploit a parametrized form |22j of the Ar- 
gonnc Vl8 + 5v + UIX* model of the EoS given in 
|23| , which is based on the most recent models for 
the nucleon-nucleon interaction with the inclusion 
of a parameterized three-body force and relativis- 
tic boost corrections. All details concerning the 
EoS and the hadronic cooling processes are given 
in . As we have seen in the previous section, it 
may be sufficient to discuss only two-flavor quark 
matter for applications to compact stars. We will 
focus on the model of the quark EoS developed 



Figure 4. Mass - radius relations for compact star 
configurations with different EoS: purely hadronic 
star with HHJ EoS (dashed), stable hybrid stars 
for HHJ - INCQMcauss EoS with 2SC (solid) and 
without 2SC phase (dash-dotted). For compar- 
ison, observational constraints on the compact- 
ness are given from the " small" compact star RX J 
1856 and from the high surface redshift object 
EXO 0748 which can both be obeyed by our hy- 
brid star EoS. 

In Refs. [221 it has been demonstrated on the 
example of the hadron-quark mixed phase that 
finite size effects might play a crucial role by sub- 
stantially narrowing the region of the mixed phase 
or even preventing its appearance. Therefore we 
omit the possibility of the hadron-quark mixed 
phase in our model where the quark phase arises 
by the Maxwell construction. For the case of 
the instantaneous Gaussian formfactor with 2SC 
phase the quark core appears for M > 1.21 M Q . 
Without 2SC phase or for the Lorentzian or NJL 
formfactor no stable hybrid stars are obtained, 
see Fig. 0] 

In order to describe the cooling of stable hybrid 
star configurations, we use the approach to the 
hadronic cooling developed recently in Ref. |2U 
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and add the contribution due to the quark core, 
see jS] . Once in the quark matter core of a hybrid 
star unpaired quarks are present the dominant 
quark direct Urea (QDU) process is operative and 
cools the star too fast in disagreement with the 
observational data ■ If one allows for a residual 
weak pairing channel as, e.g., the CSL one with 
typical gaps of A - 10 keV - 10 MeV, see j2Z], 
the hybrid star configuration will not be affected 
but the QDU cooling process will be efficiently 
suppressed as is required for a successful descrip- 
tion of compact star cooling. Since we don't know 
yet the exact pairing pattern for which also the 
residual (in 2SC unpaired) quarks are paired in 
the 2-flavor quark matter, we will call this hypo- 
thetical phase "2SC+X" . In such a way all the 
quarks may get paired, either strongly in the 2SC 
channel or weakly in the X channel. If now the 
X-gap is of the order of 1 MeV, then the QDU 
process will be effectively switched off and the 
cooling becomes too slow, again in disagreement 
with the data This seems to be a weak point 
and we would like to explore whether a density- 
dependent X-gap could allow a description of the 
cooling data within a larger interval of compact 
star masses. For such a density-dependent X-gap 
function we employ the ansatz 

Ax(m) = Ac exp[-a(^ - /l« c )/Mc] . (10) 

In Fig. [S]we show the resulting cooling curves for 
the parameters A c = 1.0 MeV and a = 10. We 
observe that the mass interval for compact stars 
which obey the cooling data constraint ranges 
now from M — 1.32 M & for slow coolers up to 
AI = 1.75 Mq for fast coolers such as Vela. 

5. Conclusion 

We have presented the QCD phase diagram 
obtained within the NCQM approach. We find 
that in the framework of our model for neutron 
star applications it is sufficient to consider only 
two quark flavors since the critical densities of 
strange quark deconfinement are not accessible 
in typical compact star interiors. We investi- 
gate the phase transition between superconduct- 
ing quark matter and hadronic matter under neu- 
tron star constraints and find that stable hybrid 
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Figure 5. Cooling of hybrid star configurations of 
Fig. 0]with color superconducting quark matter 
core in 2SC+X phase. Different lines correspond 
to hybrid star masses in units of the solar mass. 

stars are possible with large quark matter cores 
in the 2SC phase. The solution of the cooling 
evolution shows that the presence of ungapped 
quark species entails too fast cooling due to the 
QDU process in disagreement with modern neu- 
tron star cooling data. A successful description 
of hybrid star cooling requires a pairing pattern 
in the quark matter core (if it exists) where all 
quarks are gapped and a weak pairing channel is 
present with a gap which does not exceed about 
10 — 50 keV with a decreasing density dependence. 
A precise microscopic model for this conjectured 
"2SC+X" phase is still lacking. Possible candi- 
date for a weak pairing channel is the color- spin- 
locking (CSL) phase. We have demonstrated on 
the example of hybrid star cooling how astrophys- 
ical observations could provide constraints for mi- 
croscopic approaches to the QCD phase diagram 
in the domain of cold dense matter not accessi- 
ble to heavy-ion collision experiments and lattice 
QCD simulations. 
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